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Abstract 

Our purpose is to present all static solutions of the Goldstone model 
£T) | on a circle in 1 + 1 dimensions with an antiperiodicity condition imposed on 

the scalar fields. Jacobi elliptic and standard trigonometric functions are 
used to express the solutions found and stability analysis of the latter is 
what follows. Classically stable quasi-topological solitons are identified. 
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1 Introduction 

The purpose of this note is to present the complete list of static classical solu- 
Qh! tions of the Goldstone model on a circle S 1 of radius L and the corresponding 

bifurcation tree together with a study of the stability of our solutions. It comes 
as a supplement to a previous note [1] which was searching for such solutions 
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but with periodicity condition imposed there. Many of the results have obvious 
similarities with those in [1]. We notice those and make comparisons with this 
reference. We also check the case of mixed boundary conditions and place a note 
with our conclusions in the end. The following analysis has both mathematical 
[8] and physical interest as it can be useful for the search of stable solitons in 
the two-Higgs standard model (2HSM) or the minimal supersymmetric standard 
model. 



2 The classical solutions 

The Langrangian of our model is 

^ = ^(^l) 2 + ^(^2) 2 -U(0 1 ,0 2 ), /I = 0,1 (1) 

with 0i and (j>2 real Higgs fields and potential of the following form: 

U(0!,0 2 ) = i(0? + ^-l) 2 (2) 
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We impose antiperiodic boundary conditions on the scalar fields with coordinates 
x E [0,2nL]. 

<Pi (x + 2ttL) = -0, (x) , i = 1, 2 (3) 

The general static solution of model (1) can be expressed in terms of Jacobi 
elliptic functions: 

0i = VR cos Q, 02 = V^R sin Q (4) 

= ai + a 2 sn 2 (v / 2 A(rr - x ) , fc) (5) 




where C (ai (A;, A) , a 2 (A;, A)) ,x ,£ are constants while sn denotes the Jacobi el- 
liptic function; sn(z,k), sn 2 (z,k) are periodic functions on the real axis with 
periods AK(k) and 2K(k), respectively. K{k) is the complete elliptic integral 
of the first kind. Inserting (4)-(6) into the equations corresponding to (1) leads 
to several conditions on the parameters and as a consequence to three different 
types of non-trivial solutions. We start with the simplest, the trivial solution. 

2.1 The trivial solution 

The energy functional for static configurations is given by 

Apart from the vacuum solutions which have E vac = 0, one can immediately 
think of the simplest solution, the trivial one 

Ln 

0x = 2 = , E = — (8) 

which exists for all values of L. The corresponding small oscillation eigenmodes, 
labeled by j, have 

^(j) = ±aj + l/2) 2 -L 2 ), j = 0,1, 2... (9) 

This solution is stable until L = 1/2 because oo 2 (0) < for L > 1/2. Many 
additional solutions, some of which were discussed in [2,3] bifurcate from the 
solution 0i = 02 = at critical values of L. 
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2.2 Three types of non-trivial solutions 

Now, we present without many details the three non-trivial solutions and mention 
their similarities with those in [1]. The simplest case one can think, is to set a 2 = 
so that from (5) we obtain R(x) = a± and the solution becomes 



• = \l 1 " ,/ ( " T^— I (10) 



02 = \l 1 - ( A ^ j (11) 

where iV is an integer. This solution is called type-I. One can observe that the 
above solution reduces to (8) when L = N+l/2 i.e. when one of the Cj of equation 
(9) crosses zero. The Higgs field winds (2N + l)/2 times around the top of the 
Mexican hat. It's energy is given by 

7T / _ 1 \ 7T / _ 1 



BdL,N) = z \S + -) -^{*+- 2 ) ■ (12) 

If we denote the above result which corresponds to the antisymmetric case we 
study here with Ef and the result of [1] for the same type of solutions with Ef 
where S means "Symmetric" and A "Antisymmetric" then we notice that 

Ef(N,L) = Ef(N+±,L 

Another choice is to set a\ = on (5) so the solution now involves the Jacobi 
elliptic functions. C becomes zero as well and the solution is 

0! = 2kAsn (y/2Ax,kj , 2 = , A 2 = ^ l - — (13) 

This solution is called type-II. In fact, it corresponds to an oscillation of the Higgs 
field in the 2 = plane about the origin 0i = = 2 - If we take account of the 
antiperiodicity condition (3), then the argument k of the Jacobi elliptic function 
is related to the radius L of S 1 through the following formula 
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Figure 1. The solutions of equations (21) and (22) are plotted as functions of k for two values of (n,m). 
The dashed curve indicates the limit (20). 

L= ( 2m+1)Jf(t)v _ 

7T 

for some integer m. When we reach the limit of k — > (i.e. L — > m + |) the 
solution (13) approaches (8). The energy of the solution above is given by the 
integral 



E 



n 



4 (2m + 1) 
v^A(l + fc 2 ) : 



K(k) 



dy 



k sn (y, k) — 



l + k 2 \ 2 . 2k 2 -l-k 4 ' 



+ 



and by means of Elliptic integrals it becomes 



E 



n 



4(2m + 1) 
v/2A (1 + k 2 f 



K{k) (3k 4 + 2k 2 -5) + ^ (k 2 + 1) 
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Two specific values of k with the corresponding results follow 
for k = 



(15) 



E u (L = m+l/2,m) 



(2m + 1) 7T 



(16) 
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for k = 1 



B „(L = 00 , m) = l^ + i) (17) 

The comparison of E 1 ^ which represents our solutions, with Ef T which represents 
the symmetric case of the same type of solutions studied in [1], implies that 

Etj{m) = EUm + \) 

If none of a±, a 2 are zero then we are led to type-Ill solutions where we have 
the following conditions for a±, a 2 and C 2 

ai = ^(l-2A 2 (l + £; 2 )), a 2 = 4k 2 A 2 (18) 
3 

C 2 =i- (1 + (4A; 2 - 2) A 2 ) (l + (4 - 2A; 2 ) A 2 ) (l - 2A 2 (l + A; 2 )) 

=1 (1 + (4k 2 - 2) A 2 ) (1 + (4 - 2A; 2 ) A 2 ) a, (19) 
9 

Here, one can explicitly observe the fact that when a± — then C 2 = as 
well. In addition, R and C 2 must not be negative. This means another condition 

A 2 < (20) 

" 2 (1 + k 2 ) 1 ' 

with the equality leading to type-II solutions. 

In order for Q to satisfy Q (x + 2nL) = Q (x) + it and R to be periodic (so for 
the solution to be antiperiodic as we want) on [0, 2nL] the following equations 
must hold: 

c l W) dy = (2n+1) * (21) 

^ _ (2m + 1) Kjk) 

v^vrA V ' 

where m, n are positive integers which respectively determine the number of os- 
cillations of the modulus of the Higgs field and the number of the Higgs field 
windings around the origin 0i = = 02 in a period 2ixL. 

Now we return to (18) and (19) and analyze further type-Ill solutions. Solving 
these equations for k = 0, remembering that K(0) = ir/2, we find the critical 
values of L where these solutions start to exist: 



A 2 = ™2 - - =► L 2 = ^(2n+1) 2 -^ (23) 



J= (2n+l) 2 
6(2n+l) 2 -4m 2 4 V ' 2 
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The expression for A 2 together with (20) leads to the condition 2n + 1 > m on 
the integers m and n (m and n^O) . For any n and m ^ 2n + 1 the coefficient 

s\V cn^ in 4" 1~» si i n+ *"irrr*o 1 




Figure 2. The energies of some of the solutions are plotted as functions of the parameter L. The four stars 
on the upper branch indicate the four bifurcation points on the N=2 type-I solution. The two lower stars 
indicate the two bifurcation points on the N=l type-I solution. The two numbers in parentheses refer to the 
number of negative and zero modes of the corresponding solution. 

(6) is proportional to k 2 and if one expands the solution in powers of k 2 is led 
to the following formulae 



A 2 



m 



6 (2ra + 1)' - 4m 2 



1 + 



k' : 



<D(k A ) 



L 2 = \{2n + l) 2 ^ + 0{k 4 ) 



(24) 
(25) 



C 2 = 



4 ((2ra + l) 2 -m' 



n 



(3 (2t 



1^ 



2m 2 )* 



+ 0(k 4 ). 



(26) 



This expansion is necessary as there is no closed form for integral (21). In the 
limit k = we observe that the 2n solutions of type-Ill yield the type-I solution 
with N = n + 1/2 at L 2 = (3/4) (2n + l) 2 - m 2 /2 , where m = 1, 2, 2n. For 
fixed values of k, n, m we find a single solution A 2 (k, (2n + 1) /m) obeying the 
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following properties 



A 2 (k = 0, (2n + 1) /m) = — A 2 (fc = 1, (2n + 1) /m) = 1/4 

6 (2n + 1) — 4m 2 

(27) 

This is illustrated in figure 1 for two different values of the pair (n, m) by the 
solid curves, the dashed curve representing the limit (20). 

The energy of the general type-Ill solution is given by the integral 



E 



in 



2m + 1 

K(k) 



I 



dy 



( Ql - 1 + a 2 sn 2 (y, k)f + ± (l + 16A 4 (k 2 - 1 - A; 4 )) 



(28) 



which, by means of K(k), E{k) becomes 



E 



2m + 1 



in 



y/2 A 

For A; = 1 one obtains 



K(k) ( - A 2 k 2 (1 + A 2 - A 2 k 2 ) - - A 2 (2 + A 2 ) + —) + - A 2 E(k) 
\9 9 18/ 3 



4 (2m + 1) 
Em (1, m, n) = ; = E n (1, m) . 



(29) 



This shows that solution III approaches solution II in the limit k — > 1, a fact 
which also appears in [1] . 

The energies of the solutions of some low-lying branches are presented in figure 
2 as functions of L . For L > 5/2 all four types of solutions coexist and satisfy 



E I (L, N = 0) < Ei {L, N = 1) < E m (L, m = 1, n = 1) < £// (L, m = 1) < E (L) . 

(30) 

The four stars on the upper part of the figure show the position of the four 
bifurcation points L = 3.279, 3.775,4.093,4.272 of the n = 2, m = 3,2,1,0 
solutions respectively from the N = 2 type-I solution. The two lower stars 
show the bifurcation values (L = 2.179, 2.5) of the n = 1, m = 2, 1 solutions 
respectively from the N — 1 type-I solution. 



3 Stability Analysis 
3.1 Type-I solutions 

In order to analyze the stability of type-I solutions we follow the steps done in 
[1] which uses notions that can be found in [10] and [15]. Thus perturbing the 
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fields (f) a ,a — 1,2 around the classical solution (10)-(11), denoted here by 4>t , 

<Pa (X) = <Pa (X) + 7] a (x) exp (-iut) (31) 

we are led to the following equation for the normal modes: 

(32) 





a(n,l) e 

where to 2 is the eigenvalue and c = cos (^Nx/lJ , s = sin (^Nx/Lj , iV+1/2 ee JV. 

The complete list of eigenvalues of the operator A (^N, Lj for N > 1 can be 
obtained by classifying its invariant subspaces. This can be done by using the 
Fourier decomposition. Type-I solution has a twisted translational invariance and 
one can check that for an integer n > N the following finite-dimensional vector 

spaces 1 are preserved by A ^N, Lj 

V n = Span |a p cos — , (3 P sin — , p — n = \ mod2Nj , \p\ < nj , (34) 



~ ( px px ( ~\ ~\ 

V n = Span |a p cos — , — f3 p sin — , p — n = ymod2Nj , \p\ < nj (35) 

under the following condition 

a p = P P if p - n + 2N > (36) 

where a p , (3 P are arbitrary constants. The operator A (jV, can then be diag- 

onalized on each of the finite-dimensional vector spaces above, leading to a set of 
algebraic equations. 

To be more specific, define Ai = 2 (1 — A2), A2 = (jV/Lj . Also, consider the 



vector 



1 These vector spaces should have been the same as in (34) of [1] (with N — > N) but they 
are not, due to a misprint in [1]. Here we correct this by writing these spaces explicitly in 
equations (34) and (35). 
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Acting with the operator A (n, on one of the above vectors (say V^ +k ) we 
have the following steps: 



Ai / 2c 2 2sc 



S S 0- 



dx 2 2 V 2sc 2s 2 




Figure 3. The values (9) are plotted as functions of L for j=0,l,2 (solid curves), together with the values 
00^(1,1,-1), co 2 2 (l,0 -1), co 3 2 (l,0,l), co 4 2 (l,2 -1), co 5 2 (l,l,l), a> g 2 (l,3,-l) of (37) (dashed and dotted curves). 
The numbers indicate the multiplicity of the eigenvalues. 



N+k 



0. 



Where 1 is the 2x2 unit matrix and 



M 




sin 

-cos(^ 



which has the action 



The same happens if we choose Vf ! _ k instead. In that case, the equation above 
is MVfr_ k = Vfj +k as expected. A matrix which has exactly the same action as 

above can be found easily and this is ( ^ ^ 1 ■ Finally, we observe that 
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which implies the condition that the determinant of the above matrix must be 
zero in order to acquire the non-zero eigenvalues. 
The eigenvalues of A ( N, L ) on 14 are 



u 2 (N,k,±l) 



L 2 



L — [ N + 



1 



± 



L 2 



L 2 — 




+ 4k 2 [N + 



(37) 



for k = 0,1, 2, ... and similarly on V n for A; = 1,2,3,... . Notice that if we replace 
N + 1/2 by N then the above result is the same with eq.(36) of [1]. 

For N having a specific value and L slightly greater than N , the solution 
(10)-(11) possess 4iV negative modes corresponding to k — 1,2,3, ...,2N in (37) 
(minus sign). When L increases, the number of positive eigenmodes increases as 
well 



l^N+lf- — 
4 V ; 2 



m 



1,2, 



,2N, 



(38) 



i.e. (cf(23)) at those values of L where the solutions of type-Ill with n = N 
bifurcate from the solution of type-I. For 



L 2 > L\ r (iV) ee - (2N + If 



(39) 



1 

4 v-- ■ -/ 2 

all the modes are positive and (10)-(11) are classically stable solitons. These re- 
sults are illustrated in figure 2 for N — 1 and N = 2. The numbers in parentheses 
represent the number of negative and zero modes of the corresponding branch. 



3.2 Type-II and Type-Ill solutions 

The equations for the fluctuations about the solution (13) decouple to take the 
form of Lame equations: 



dy 2 

d 2 
dy 



+ Qk 2 sn 2 (y, k) j r/i = Vt\r]i 
2 + 2k 2 sn 2 (y, k) \ 772 = fi 2 r/ 2 



(40) 
(41) 
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where 



x = VTT&y, Q 2 a =(u 2 a + l)(k 2 + l), a = 1,2 (42) 

and u 2 is the effective eigenvalue of the relevant operator. Equations (40) and (41) 
admit two and one algebraic modes, respectively, with corresponding eigenvalues 

Vi\ : 4 + k 2 , 1 + Ak 2 and Q 2 2 : 1 + k 2 (43) 

The corresponding values of a; 2 follow immediately from (42); they have signature 
(+,+) and (0), respectively. 

It's a property of the Lame equation that the solutions determined alge- 
braically correspond to the solutions of the lowest eigenvalues. The remaining 
part of the spectrum therefore consists of positive eigenvalues. The spectrum of 
the equation (40) was studied perturbatively in [9], while the relation between 
the Lame equation and the Manton-Samols sphalerons was first pointed out in 
[11-13]. 

Concerning the stability equation of type-Ill solution, it seems to be more 
difficult to deal with as the presence of the Jacobi Elliptic functions prevent us 
from having an analytic expression for Vl{x) in (4). Thus, we are not able to 
follow the steps done in the case of type-I solutions (i.e. find the M matrix) 
where trigonometric functions are easier to handle. Detailed analysis especially 
for the case of Lame equations can be found in [12] and [15]. 

4 Conclusion 

A detailed analytical study of the static solutions of the Goldstone model on a 
circle has been given in [1] and we followed the same path here for our boundary 
conditions. Many results of [1] are connected with ours by a simple change 
on variables used. We write them explicitly above whenever is necessary. We 
also note our effort to impose mixed boundary conditions as well. Specifically, 
we enforced 0i to be antiperiodic and 02 periodic but there was no solution to 
satisfy this choice so it's needless to extend beyond this small remark here. 

Many details were presented on the stability analysis on section 3 for the 
solutions found. Classically stable solitons were identified, together with the 
range of the parameter L for which they are stable. 

This simpler model we present here and many others can be connected and 
can also give us the experience to deal with realistic (3+l)-dimensional particle 
physics models in our search for possible metastable localized solitons. Ref. [1] 
on which this note was based, has connections with [9] and [10] as it concerns 
the branches of their solutions, also with [2-6], while there are also interesting 
physical applications [14] as it is already mentioned in [1]. Soliton solutions 
in (3+l)-dimensional models with antiperiodic boundary condition imposed on 
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one spatial dimension, which is compactified on S , is analyzed in [16-18] where 
supersymmetry breaking is found. 
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Our purpose is to present all static solutions of the Goldstone model on a circle 
in 1+1 dimensions with an antiperiodicity condition imposed on the scalar helds. 
Jacobi elliptic and standard trigonometric functions are used to express the solutions 
found and stability analysis of the latter is what follows. Classically stable quasi- 
topological solitons are identified. 

1 Introduction 

The purpose of this note is to present the complete list of static classical solutions of the 
Goldstone model on a circle S 1 of radius L and the corresponding bifurcation tree together 
with a study of the stability of our solutions. It comes as a supplement to a previous 
note [1] which was searching for such solutions but with periodicity condition imposed 
there. Many of the results have obvious similarities with those in [1]. We notice those 
and make comparisons with this reference. We also check the case of mixed boundary 
conditions and place a note with our conclusions in the end. The following analysis 
has both mathematical [8] and physical interest as it can be useful for the search of 
stable solitons in the two-Higgs standard model (2HSM) or the minimal supersymmetric 
standard model. 

2 The classical solutions 

The Langrangian of our model is 



We impose antiperiodic boundary conditions on the scalar helds with coordinates x £ 



Abstract 





with (pi and <f> 2 real Higgs helds and potential of the following form: 





[0,2wL]. 



1 



fa(x + 2wL) = -<f>i(x) , i = 1,2 (3) 

The general static solution of model (1) can be expressed in terms of Jacobi elliptic 
functions: 



v^Rcosfi, ^ 2 = v / ftsinfi (4) 



R(x) = a,\ + a 2 sn 2 (\/2 A(x — x ) , k) (5) 



n(l) = °h Wf* (6) 

where C (ai (A;, A) , a 2 A)) ,x ,^ are constants while sn denotes the Jacobi elliptic func- 
tion; sn(z,k), sn 2 (z,k) are periodic functions on the real axis with periods 4K(k) and 
2K(k), respectively. K{k) is the complete elliptic integral of the first kind. Inserting (4)- 
(6) into the equations corresponding to (1) leads to several conditions on the parameters 
and as a consequence to three different types of non-trivial solutions. We start with the 
simplest, the trivial solution. 

2.1 The trivial solution 

The energy functional for static configurations is given by 

Apart from the vacuum solutions which have E vac = 0, one can immediately think of the 
simplest solution, the trivial one 



fa = fa = , E = — (8) 

which exists for all values of L. The corresponding small oscillation eigenmodes, labeled 
by j, have 

^ 2 (j) = ^((j + l/2) 2 -£ 2 ), J = 0,1,2... (9) 

This solution is stable until L = 1/2 because cu 2 (0) < for L > 1/2. Many additional 
solutions, some of which were discussed in [2,3] bifurcate from the solution fa = <f> 2 = 
at critical values of L. 
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2.2 Three types of non-trivial solutions 

Now, we present without many details the three non-trivial solutions and mention their 
similarities with those in [1]. The simplest case one can think, is to set a 2 = so that 
from (5) we obtain R(x) = a\ and the solution becomes 



where N is an integer. This solution is called type-I. One can observe that the above 
solution reduces to (8) when L = N + 1/2 i.e. when one of the u of equation (9) crosses 
zero. The Higgs field winds (2iV + l)/2 times around the top of the Mexican hat. It's 
energy is given by 



E , {L , N) = l( N + ^ 2 -JL( N + ^\ (12) 

If we denote the above result which corresponds to the antisymmetric case we study here 
with Ef and the result of [1] for the same type of solutions with Ef where S means 
"Symmetric" and A "Antisymmetric" then we notice that 

Ef{N,L) = Ef(N + ±,L) 

Another choice is to set a\ = on (5) so the solution now involves the Jacobi elliptic 
functions. C becomes zero as well and the solution is 



^ = 2k\sn (V2Kx^ , <^2 = 0, A 2 = 2 ^\ k2) (13) 

This solution is called type-II. In fact, it corresponds to an oscillation of the Higgs field in 
the ^2 = plane about the origin cf>i = = <f> 2 - If we take account of the antiperiodicity 
condition (3), then the argument k of the Jacobi elliptic function is related to the radius 
L of S 1 through the following formula 
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Figure 1. The solutions of equations (21) and (22) are plotted as functions of k for two values of (n,m). 
The dashed curve indicates the limit (20). 



z = (2m + 1)A - Wy _ 

7T 



(14) 



for some integer m. When we reach the limit of k — > (i.e. L — > m + |) the solution 
(13) approaches (8). The energy of the solution above is given by the integral 



E 



4 (2m + 1) 



ii 



K(k) 



dy 



k sn (j/, k) 



1 + k 



2\ 2 



V2A(l + F)Vo 
and by means of Elliptic integrals it becomes 



+ 



2k 2 - 1 - k 4 



(15) 



E 



ii 



4 (2m + 1) 
V2A(1 + F 



K{k) (3F + 2k 2 - 5) + ^ (* 2 + i; 
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Two specific values of k with the corresponding results follow 
for k = 



= m + 1/2, m) 



(2m + 1) 7T 



for k = 1 



(16) 



£7/ (1/ = 00, m) 



4 (2m + 1) 
3^2 



(17) 
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The comparison of Efj which represents our solutions, with Efj which represents the 
symmetric case of the same type of solutions studied in [1], implies that 

Efjim) = E s n (m+ 1 -) 

If none of ai,a 2 are zero then we are led to type-Ill solutions where we have the 
following conditions for ai,a 2 and C 2 

ai = - (l — 2A 2 (l + k 2 )) , a 2 = 4k 2 A 2 (18) 



C 2 =_L (l + (4k 2 - 2) A 2 ) (1 + (4 - 2k 2 ) A 2 ) (1 - 2A 2 (l + k 2 )) 

= 1 (i + ( 4 £2 _ 2 ) a 2 ) (1 + (4 - 2k 2 ) A 2 ) ai (19) 

Here, one can explicitly observe the fact that when a\ = then C 2 = as well. In 
addition, R and C 2 must not be negative. This means another condition 



A 2 < ptt (20) 

~ 2(1 + k 2 ) K ' 

with the equality leading to type-II solutions. 

In order for O to satisfy O (x + 2ttL) = O (x) + tt and R to be periodic (so for the 
solution to be antiperiodic as we want) on [0, 2ttL] the following equations must hold: 

C T L ^dy = (2n + l)7r (21) 

Jo R{y) 



(2m + 1) K(k) , , 

L = ^J- — — 22 

where m, n are positive integers which respectively determine the number of oscillations 
of the modulus of the Higgs field and the number of the Higgs field windings around the 
origin cf>i = = <f> 2 in a period 2ttL. 

Now we return to (18) and (19) and analyze further type-Ill solutions. Solving these 
equations for k = 0, remembering that K(0) = 7r/2, we find the critical values of L where 
these solutions start to exist: 



A 2 = 5 L 2 = -(2n + l) 2 23 

6(2n + l) 2 -4m 2 4 V ' 2 K ' 

The expression for A 2 together with (20) leads to the condition 2n + l > m on the integers 
m and n (m and n ^ 0) . For any n and m ^ 2n + 1 the coefficient of sn 2 in the integral 
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L 

Figure 2. The energies of some of the solutions are plotted as functions of the parameter L. The four stars 
on the upper branch indicate the four bifurcation points on the N=2 type-I solution. The two lower stars 
indicate the two bifurcation points on the N=l type-I solution. The two numbers in parentheses refer to the 
number of negative and zero modes of the corresponding solution. 

(6) is proportional to k 2 and if one expands the solution in powers of k 2 is led to the 
following formulae 



A 



m 



6(2n + 1) - 4m ; 



1 + ^ )+o(k*) 



(24) 



L 2 = 3 -(2n + lf-^- + 0(k*) 



(25) 



, 4 ((2n + l) 2 - m 2 ) 2 (2n + 1? ... 
3 ( 3 (2n + l) 2 -2m2) 3 V ' 

This expansion is necessary as there is no closed form for integral (21). In the limit k = 
we observe that the 2n solutions of type-Ill yield the type-I solution with N = n + 1/2 
at L 2 = (3/4) (2n + 1) - m 2 /2 , where m = 1,2, ...,2n. For fixed values of k, n, m we 
hnd a single solution A 2 (k } (2n + 1) /m) obeying the following properties 
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m 



A 2 (k = 0, (2n + 1) /m) = — — — ^ — , A (k = 1, (2n + 1) /m) = 1/4 (27) 



6(2n + 1) - 4m 

This is illustrated in figure 1 for two different values of the pair (n, m) by the solid curves, 
the dashed curve representing the limit (20). 
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The energy of the general type-Ill solution is given by the integral 



E 



in 



2m + 1 

K(k) 



dy 



( ai - 1 + a 2 sn 2 (y } k)) + - (l + 16A 4 (k 2 - 1 - k 4 )) 



(28) 



which, by means of K(k) } E{k) becomes 



E 



2m + 1 



in 



V2~A 

For k = 1 one obtains 



K(k) l-A 2 k 2 (l + A 2 -A 2 k 2 " 



-A 2 (2 + A 2 ) + — 
9 v > 18 



+ -A 2 £(£) 



E m (1, m, n) = 4(2 3 ^ 1} = (1, ™) • (29) 

This shows that solution III approaches solution II in the limit k — > 1, a fact which also 
appears in [1] . 

The energies of the solutions of some low-lying branches are presented in figure 2 as 
functions of L . For L > 5/2 all four types of solutions coexist and satisfy 



Ei (L, N = 0) < E I (L, N = 1) < Em (L, m = 1, n = 1) < E u (L, m = 1) < E (L) . 

(30) 

The four stars on the upper part of the figure show the position of the four bifurca- 
tion points L = 3.279, 3.775,4.093,4.272 of the n = 2, m = 3, 2, 1, solutions respec- 
tively from the N = 2 type-I solution. The two lower stars show the bifurcation values 
(L = 2.179, 2.5) of the n = 1, m = 2, 1 solutions respectively from the N = 1 type-I 
solution. 



3 Stability Analysis 
3.1 Type-I solutions 

In order to analyze the stability of type-I solutions we follow the steps done in [1] which 
uses notions that can be found in [10] and [15]. Thus perturbing the fields <f> a , a = 1,2 
around the classical solution (10)-(11), denoted here by (f) c J , 



<f> a {x) = <Pa (x) + r] a (x)exp(-iut) (31) 
we are led to the following equation for the normal modes: 

Vi\ 2 ( Vi 



A [N,L)[ - I= U M - ) (32) 
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A(N,L)=-— + 2|l-[-] | ( 1 - )-(^\ (33) 



where to 2 is the eigenvalue and c = cos ( Nx/ L ) , s = sin ( Nx/ L ) , N + 1/2 = iV. 





The complete list of eigenvalues of the operator A (^iV, for N > 1 can be obtained 
by classifying its invariant subspaces. This can be done by using the Fourier decomposi- 
tion. Type-I solution has a twisted translational invariance and one can check that for an 
integer n > N the following finite-dimensional vector spaces 1 are preserved by A ( N, L 



( px px / ~ \ *| 

V n = Span |a p cos — , f3 p sin — , p — n = ymod2Nj , \p\ < nj , (34) 

~ ( px px / - \ 'I 

V n = Span |a p cos — , —f3 p sin — , p — n = ymod2N J , \p\ < n j (35) 

under the following condition 

a p = p p if p - n + 2N > (36) 



where a v , f3 p are arbitrary constants. The operator A (^iV, can then be diagonalized on 
each of the finite-dimensional vector spaces above, leading to a set of algebraic equations. 
To be more specific, define Ai = 2 (1 — A 2 ), A 2 = ( N/L ) . Also, consider the vector 



V*- t ) W+ * - \ S in (S±A x J ■ W -* - [ S in * 



Acting with the operator A yN } Lj on one of the above vectors (say V^ +k ) we have the 
following steps: 

d 2 Ai ( 2c 2 2sc \ 2 



dx~ 2+ Y 2sc 2s 2 ~ X2 ~" \ V N +k 



1 These vector spaces should have been the same as in (34) of [1] (with N —> N) but they are not, due 
to a misprint in [1]. Here we correct this by writing these spaces explicitly in equations (34) and (35). 
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Figure 3. The values (9) are plotted as functions of L for j=0,l,2 (solid curves), together with the values 
00^(1,1,-1), co 2 2 (l,0 -1), co 3 2 (l,0,l), co 4 2 (l,2 -1), co 5 2 (l,l,l), co 6 2 (l,3,-l) of (37) (dashed and dotted curves). 
The numbers indicate the multiplicity of the eigenvalues. 



d 2 Ax 



' N+k 



0. 



Where 1 is the 2x2 unit matrix and 



M 



cos (^) 
sin 



sm 



cos 



2Nx 
L 

( 2Nx 
\ L 



which has the action 



The same happens if we choose Vj^_ k instead. In that case, the equation above is 
MV)y_ fc = Vfi +k as expected. A matrix which has exactly the same action as above 

can be found easily and this is (^ ^ ^ ^ • Fi na Uy? we observe that 

0¥) 2 +^-a 2 -- 2 ^ 

Ai fN-kV i Ai ^ ,2 I V *+ k 



1L 

2 \ L I 1 2 



+ ^ - A 2 - u> 



which implies the condition that the determinant of the above matrix must be zero in 
order to acquire the non-zero eigenvalues. 
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The eigenvalues of A I N, L J on V n are 



u 2 (N,k,±l)= — 



L 2 



N + -) +k 2 



± 1 2 \ 



L 2 




+ 4k 2 [ N + - 



(37) 



for k = 0, 1, 2, ... and similarly on V n for A; = 1,2,3,... . Notice that if we replace N + 1/2 
by N then the above result is the same with eq.(36) of [1]. 

For N having a specific value and L slightly greater than iV, the solution (10)-(11) 
possess 4iV negative modes corresponding to k = 1, 2, 3, 2iV in (37) (minus sign). 
When L increases, the number of positive eigenmodes increases as well 



r -> 3 . , T , 2 m 
L 2 = ^(2N + 1) 2 -—, 



m = 1,2, ...,27V, 



(38) 



i.e. (cf(23)) at those values of L where the solutions of type-Ill with n = N bifurcate 
from the solution of type-I. For 



L 2 > L 2 cr (N) = - (2JV + 1) 



(39) 



all the modes are positive and (10)-(11) are classically stable solitons. These results are 
illustrated in figure 2 for N = 1 and N = 2. The numbers in parentheses represent the 
number of negative and zero modes of the corresponding branch. 

3.2 Type-II and Type-Ill solutions 

The equations for the fluctuations about the solution (13) decouple to take the form of 
Lame equations: 



1 
dy 



2 + 6k 2 sn 2 (y, k) \ r] x = fi^i 



(40) 



where 



dy 



2 +2k 2 sn 2 (y,k) \r, 2 = 0^ 2 



(41) 



Vl + k 2 



y 



ill EE (u 2 a + 1) (k 2 + 1) 



1,2 



(42) 



and u 2 a is the effective eigenvalue of the relevant operator. Equations (40) and (41) admit 
two and one algebraic modes, respectively, with corresponding eigenvalues 
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ttl : 4 + k 2 , 1 + 4k 2 



and 



tt 2 2 : 1 + k 2 



(43) 



The corresponding values of u 2 follow immediately from (42); they have signature (+,+) 
and (0), respectively. 

It's a property of the Lame equation that the solutions determined algebraically cor- 
respond to the solutions of the lowest eigenvalues. The remaining part of the spectrum 
therefore consists of positive eigenvalues. The spectrum of the equation (40) was studied 
perturbatively in [9], while the relation between the Lame equation and the Manton- 
Samols sphalerons was first pointed out in [11-13]. 

Concerning the stability equation of type-Ill solution, it seems to be more difficult 
to deal with as the presence of the Jacobi Elliptic functions prevent us from having an 
analytic expression for Q(x) in (4). Thus, we are not able to follow the steps done in the 
case of type-I solutions (i.e. find the M matrix) where trigonometric functions are easier 
to handle. Detailed analysis especially for the case of Lame equations can be found in 
[12] and [15]. 

4 Conclusion 

A detailed analytical study of the static solutions of the Goldstone model on a circle 
has been given in [1] and we followed the same path here for our boundary conditions. 
Many results of [1] are connected with ours by a simple change on variables used. We 
write them explicitly above whenever is necessary. We also note our effort to impose 
mixed boundary conditions as well. Specifically, we enforced cf>i to be antiperiodic and <f> 2 
periodic but there was no solution to satisfy this choice so it's needless to extend beyond 
this small remark here. 

Many details were presented on the stability analysis on section 3 for the solutions 
found. Classically stable solitons were identified, together with the range of the parameter 
L for which they are stable. 

This simpler model we present here and many others can be connected and can also 
give us the experience to deal with realistic (3+l)-dimensional particle physics models in 
our search for possible metastable localized solitons. Ref. [1] on which this note was based, 
has connections with [9] and [10] as it concerns the branches of their solutions, also with 
[2-6], while there are also interesting physical applications [14] as it is already mentioned in 
[1]. Soliton solutions in (3+l)-dimensional models with antiperiodic boundary condition 
imposed on one spatial dimension, which is compactified on S 1 , is analyzed in [16-18] 
where supersymmetry breaking is found. 
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